Temperature-dependent spin resonance energy in iron pnictides and multiband s± 

Eliashberg theory 



G.A. UmmarincQ 
Dipartimento di Fisica and C'NISM, Politecmco di Torino, 
Corso Duca degli Abruzzi 24, 10129 Torino, Italy 

The phenomenology of iron-pnictides superconductors can be explained in the framework of a 
three bands s± wave Eliashberg theory with only two free parameters plus a feedback effect i.e. the 
effect of the condensate on the antiferromagnetic spin fiuctuactions responsible of the superconduc- 
tivity in these compounds. I have examined the experimental data of four materials: LaFeAsOi_xFx, 
SmFeAsOi_xFx, Bai_xKxFe2As2, and Ba(FexCoi_x)2As2 and I have found that it is possible to re- 
produce the experimental critical temperature and gap values in a moderate strong-coupling regime: 
\tot ~ 1.7-2.0. 

PACS numbers: 74.70.Dd, 74.20.Fg, 74.20. Mn 



The new class of Fe-based compounds |T}jS] just as the 
cuprates jlj and the heavy fermions [S] have all some 
similar caracteristics. For example the high values of 
rate 2A/Tc or the presence of the pseudogap [H IHl [Z]- 
For all three class of material it is proposed the super- 
conductivity to be mediated by antiferromagnetic spin 
fiuctuactions [H [HI IH]- The most obvious difference is 
that almost all the iron compounds present a multi- 
band behavior while in HTCS and in heavy fermions 
this was detected only in some particular cases. The 
multi-band nature of Fe-based superconductors may give 
rise to a multi-gap scenario [10] that is indeed emerging 
from many different experimental data with evidence for 
rather high gap ratios, w 2 — 3 [TT]. In this regard nei- 
ther a three-band BCS model [H [TH [O] nor a four-band 
Eliashberg model [M] with small values of the coupling 
constants and large boson energies are adequate: the for- 
mer can only account for the gap ratio and but not 
for the exact experimental gap values and the latter pro- 
vides a calculated critical temperature larger than the 
experimental one. The high experimental value of the 
larger gap suggests that high values of the coupling con- 
stants might be necessary to explain the experimental 
data within a three-band model [TS] [H] : one has there- 
fore to employ the Eliashberg theory for strong coupling 
superconductors [13 IS]. In my early works [13 US] 
I found that a three-band Eliashberg model allows to 
reproduce various experimental data, this suggests that 
these compounds can represent a case of dominant nega- 
tive interband-channel superconductivity (s± wave sym- 
metry) with small typical boson energies (~ 10 meV) but 
too high values of the electron-boson coupling constants 
(1.9 < Xtot < 5.9). The way for solve this problem is 
suggested by experimental measurement of Inosov and 
coworkers [T7] : they find that the temperature evolution 
of the spin resonance energy follows the superconduct- 
ing energy gap and this should indicate a feedback effect 
[H [iHl [19] of the condensate on the spin fiuctuactions. I 
assume that this is the starting point of my argumenta- 
tion. The procedure is as follows: first of all I choose the 



experimental low temperature spin resonance as repre- 
sentative boson energy and I fix the two remaining free 
parameters to reproduce the exact experimental gap val- 
ues, then, with the same parameters, I calculate the 
critical temperature T*. I find always T* :$> T^^p where 
T^'^P is the experimental critical temperature. In the suc- 
cessive step I use the same input parameters utilized be- 
fore except for the electron-boson spectral functions that 
have an energy peak with the same temperature depen- 
dence of the superconductive gap. Of course at T = T* 
the energy peak is equal to zero while at T = K the 
new spectral functions are equal to old ones. In this way, 
taking into account the feedback effect of the condensate 
m [IHl [m on the antiferromagnetic spin fiuctuactions I 
could explain the experimental data (the gap values and 
the critical temperature) in a model which has only two 
free parameters in a moderate strong coupling regime 
{Xtot ^1-7- 2). 

I choose four representative cases (three hole type and 
one electron type): LaFeAsOi_xFx, SmFeAsOi_xFx, 
Bai_xKxFe2As2, and Ba(FexCoi_x)2As2. The electronic 
structure of the compounds hole type can be approx- 
imately described by a three-band model [^ with two 
hole bands (indicated in the following as bands 1 and 2) 
and one equivalent electron band (3) [13 [H] while for 
one electron type with one hole band (indicated in the 
following as band 1) and two equivalent electron bands 
(2 and 3) [35]. In the hole type case the s-wave order 
parameters of the hole bands Ai and A2 have opposite 
sign compared to electron band one, A3 [20] while, in 
the electron type case, Ai has opposite sign compared to 
two electron bands ones, A2 and A3 [25] . In such sys- 
tems, intraband coupling could be provided by phonons 
{ph), and interband coupling by antiferromagnetic spin 
fluctuations (sf) |20j . I summarize the experimental data 
relative to the four considered cases: 

1) the compound LaFeAsOo.gFo.i (LaFeAsOF)with 
= 28.6 K where point-contact spectroscopy measure- 
ments gave Ai(0) w 8.0 meV and A2(0) ^ 2.8 meV [3; 

2) Bao.6Ko.4Fe2As2 (BaKFeAs) with = 37 K where 
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TABLE I: The values of Qq and Xij, that allow reproducing 
the experimental gap values, are shown. Xtot is compared 
with A°ot that is the value determined in the previous works 



[151 1161 126j . In the first arrows the sf spectral functions used 
have usual shape while in the second ones have Lorentzian 
shape. 
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TABLE IL The calculated values of the gaps and of the two 
critical temperature with and without feedback effect. In the 
first arrows the sf spectral functions used have usual shape 
while in the second ones have Lorentzian shape. 



ARPES measurements gave Ai(0) = 12.1 ± 1.5 meV, 
A2(0) =5.2±1.0 meVand A3(0) = 12.8±1.4meV [21]; 

3) the compound SmFeAsOo.8Fo.2 (SmFeAsOF) with 
= 52 (Tj;^^'' = 53 K) K where, according to point- 
contact spectroscopy measurements, Ai(0) ~ 18±3 meV 
and A2(0) = 6.2 ± 0.5 meV t25j; 

4) the compound Ba(FexCoi_x)2As2 (BaFeCoAs) with 
= 22.6 K (r^''""= = 24.5 K) where, according to point- 
contact spectroscopy measurements, Ai(0) = 4.1 ± 0.4 
meV and A2(0) = 9.2 ± 1.0 meV [55]. 

is the critical temperature obtained by Andreev 
reflection measurements and T^^""= is the critical tem- 
perature extracted by transport measurements. Note 
that only in the case of ARPES the gaps are associated 
to the relevant band since point-contact spectroscopy 
measurements generally gives only two gaps, the larger 
one has been arbitrarily indicated as Ai supposing that 
Ai^lAal. 

To obtain the gaps and the critical temperature within 
the s± wave, three-band Eliashbcrg equations [22] one 
has to solve six coupled equations for the gaps Ai(iw„) 
and the renormalization functions Ziiiujn), where i is a 
band index (that ranges between 1 and 3) and uj„ are 
the Matsubara frequencies. If one neglects for simplicity 



the effect of magnetic and non-magnetic impurities, the 
imaginary-axis equations |15l I16j are: 



Zj(iw„)Aj(ia;„) = ttT V'[A,f'(iw„, ia;„) 



(2) 



where Ag(za;„, «a;,„) = Afj'(ii:j„, iw„) + A|/(ia;„, za;,„), 
A^(iw„, iuj^) = A^j{iujn, ioj^) - A^J(iw„, ia;„). 6 is the 
Heaviside function and Wc is a cutoff energy. In partic- 
ular, Afj''^^ (»..„, i^„) = 2/+°° dnna^^pp'^'^f in)/[{LOn- 
ujm)'^ + ^^]- f^iji^c) are the elements of the 
3x3 Coulomb pseudopotential matrix. Finally, 

N^{iuj,n) = Aj{iu}„i)/ ^uj^^ + A'jiiujm) and Nf{iuJm) = 
'^m/ \J<^m + A|(icj„i). The electron-boson coupling con- 



stants are defined as 



2/o 



The solution of eqs. [T] and [2] requires a huge num- 
ber of input parameters (18 functions and 9 constants); 
however, some of these parameters are related one to an- 
other, some can be extracted from experiments and some 
can be fixed by suitable approximations. As shown in 
Ref. [20], in the case of pnictides we can assume that: 
i) the total electron-phonon coupling constant is small 
[23) : ii) phonons mainly provide intraband coupling; iii) 
spin fluctuations mainly provide interband coupling. To 
account for these assumptions in the simplest way, I will 
take: Aff = Afj* = 0. (upper limit of the phonon cou- 
pling [23] « 0.35), \ll — (only interband sf coupling) 
and /i*j(wc) — Mij('^c) = [T5]. Within these approxi- 
mations, the electron-boson coupling-constant matrix Xij 
becomes: [91IT51I26]: 







Al2 Ai3 
A2I = Ai2l^l2 A23 

\ A31 = Ai3J^i3 A32 — 



(3) 



where Vij = Ni{{))/Nj{Q) and Niifi) is the normal den- 
sity of states at the Fermi level for the i-th band. In the 
hole case it is A21 — A12 = while in the electron case 
A23 = A32 = 0. In the numerical simulations I used the 
standard form for the antiferromagnetic spin fluctuaction 

[H]: a%F'p{n) = B.^n ■ n,, ■ Q{n^a. - + ^l) 

where Bij are the normalization constants necessary to 
obtain the proper values of Xij while ^ij are the peak 
energies. In all the calculations I always set Q.ij = fio- 
The maximum sf energy is ^max = lOiloi the cut-off 
energy is Uc = 30ilo and the maximum quasiparticle en- 
ergy is uJrnax = 40rio. As typical sf energy Q.q I use 
the spin resonance energy that have been measured and 
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FIG. 1: (Color online) The calculated critical temperature 
Tc with feedback effect versus standard critical temperature 
T* in three different situations: only interband sf coupling 
with standard spectral functions (black squares), interband 
sf coupling with standard spectral functions and small in- 
traband ph coupling (red circles) and only interband sf cou- 
pling with Lorentz spectral functions (dark blue triangles). 
In the bottom right insert the sf spectral function, for the 
Ba(FexCoi-x)2As2, at different temperatures (T < T^) with 
the feedback effect. 




FIG. 2: (Color online) The calculated temperature depen- 
dence of |Ai| from the solution of real axis Eliashberg equa- 
tions in the standard case (open symbol) and when the feed- 
back effect is present (solid symbol): |Ai| black squares, IA2I 
red circles and [Aaj dark blue triangles. The experimental 
data (7| i24rt26) are shown as big solid circles. 

I assume correct for all compounds examined the rela- 
tion fig = (2/5)Tc available in literature |57]. Band- 
structure calculations provide information about the fac- 
tors Vij that enter in the definition of Xij (eq. 3). In 
the case of LaFeAsOo.gFg.i I know that = 0.91 and 
1^23 = 0.53 in Bao.6Ko.4Fe2As2 h'ls = 1 and 1^23 = 2 
0, in SmFeAsOo.sF 0.2 1^13 — 0-4 and 1^23 ~ 0-5 |28j and 
in Ba(FexCoi_x)2As2 1^12 = 1-12 and 1^13 = 4.50 ^5] . 
I initially solve the imaginary-axis Eliashberg equations 



(eqs. [T]and[2]) to calculate the low-temperature value of 
the gaps (which are actually obtained by analytical con- 
tinuation to the real axis by using the technique of the 
Fade approximants) and so I fix the two free parameters 
of the model: A13 and A23 (A12). By properly select- 
ing the values of A13 and A23 (A12) it is relatively easy 
to obtain the experimental values of the gaps with rea- 
sonable values of Xtot = ^ ^n^(o)' (between 1.72 and 
2.04). However, in all the materials examined, the high 
2Ai.3/A:bT'c ratio (of the order of 8-9) makes it possi- 
ble to reproduce also the values of the large gap(s) only 
if the calculated critical temperature T* is considerably 
higher than the experimental one. For solving this prob- 
lem also present in the HTCS, I assume that exist a 
effect of feedback [U [THl [TH] of the condensate and, in 
a phenomenological way, I introduce in the Eliashberg 
equation a temperature dependence of the representative 
boson energy f)o(T) = ^^otanh{1.76^/^fJY^) that re- 
produces both the approximate gap temperature depen- 
dence in the strong coupling case [19! ^^'^ the experimen- 
tal spin resonance one ^1?) . The primary effect of this 
assumption is lowering the critical temperature leaving 
unchanged the gap values at T << T* because the crit- 
ical temperature is roughly proportional to electron bo- 
son coupling constant and to representative boson energy 
^o{T) of the material: in this case ^q{T) decreases and 
so Tc- For a completely consistent procedure it should 
used rio(r) = i}or]{T) where ri{T) is the temperature de- 
pendence part of the superfluid density p(T) = p{0)ri[T) 
and p{0) is the superfluid density at T = K. ri{T) is 
a function of Ai(iujn) and so, in this way, the numeri- 
cal solution of Eliashberg equations become remarkably 
more complex and time consuming. I am conscious that 
the temperature dependence of flo{T) is added ad hoc 
and it is not obtained self-consistent but this is an at- 
tempt in order to determine if the chosen path can lead 
to interesting results. What is important is that this 
mechanism of feedback can justify the experimental val- 
ues for the gaps, their dependence on temperature and 
the critical temperature with a model that has only two 
free parameters. Moreover, the parameters determined 
are reasonable and Xtot is very similar for all four materi- 
als examined and in agreement with the values proposed 
by other authors [21] . I solve the Eliashberg equations in 
three different situations: 1) only sf interband coupling 
is present and the sf spectral functions have usual shape; 
2) sf interband coupling with a small ph intraband con- 
tribution are present and sf spectral functions have usual 
shape; 3) only sf interband coupling is present and the sf 
spectral functions have Lorentz shape. In the first case 
the coupling constant Xtot is in the range 1.72-2.04. The 
results are almost independent from ^max because, for 
example in the case of BaFeCoAs, multiplying fimax by 
a factor two, I obtain the same values of the gaps and Tc 
with Xtot = 1-68 i.e. with a reduction of 0.18 which is 
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very small. The agreement with the experimental critical 
temperature is good. It is noticeable the small variation 
of the total coupling in the four compounds considered. 
In the second case there is also a intraband phonon con- 
tribution, equal in any band and in any compound for 
simplicity, with Af^'' = 0.35 and = 18 meV that are 
the upper limits for the ph coupling constants and the 
representative ph energies [23] . The ph spectral functions 
have Lorentzian shape [TS] with the peaks at the same 
energy: f2.y = fip'^and with half width always equal to 
2 meV { ujc = 12r2Q''). Xt^t and are practically the 
same as the previous case. This last fact indicate that 
the effect of intraband phonon contribution is negligible. 
In the third case (Lorentz shape of sf spectral functions) 
the agreement with the experimental critical tempera- 
tures is very good in all compounds but the total cou- 
pling is more large (2.38 < Xtot < 2.84). In Fig. 1 it is 
possible to see the linear relation between Tc and T* in 
all three examined cases. In table 1 are shown the inputs 
parameters of the Eliashberg equations in the first and 
third case examined for the four compounds. In tabic 
2 are shown the calculated values of the gaps and the 
critical temperatures Tc and T* obtained by numerical 
solution of Eliashberg equations. Once the values of the 
low-temperature gaps were obtained, I calculated their 
temperature dependence by directly solving the three- 
band Eliashberg equations in the real-axis formulation 
instead of using the analytical continuation to the real 
axis of the imaginary-axis solution. Of course, the re- 
sults of the two procedures are virtually identical at low 
temperature. In all cases, their behavior is rather un- 
usual and completely different from the BCS one, since 
the gaps slightly decrease with increasing temperature 
until they suddenly drop close to Tc- This arises from 
a complex non-linear dependence of the A vs. T curves 
on Xij and is possible only in a strong-coupling regime 
[53]. Curiously in all four compounds the rate T* /Tc 
is 1.27. As it is shown in Fig. 2 the calculated tem- 
perature dependencies of |Ai| are compared with the ex- 
perimental data and the agreement is very good (in the 
case of Ba(FexCoi_x)2As2 I compare the temperature de- 
pendence of the gaps with these particular experimental 
values [^: Ai = 3.8 meV and Aj = 8.2 meV and I find 
Ai2 = 0.77, Ai3 = 1.05 and = 23.43 K). In conclusion, 
I have shown that a simple Eliashberg three-band model, 
with antiferromagnetic spin fluctuactions-electrons cou- 
pling, in moderate strong-coupling regime with only two 
free parameters and a feedback effect can reproduce, in a 
quantitative way, the experimental critical temperature 
and the amplitude of the energy gaps. 

I thank R.S. Gonnelli, E. Cappelluti, Lara Benfatto 
and Sara Galasso for the useful and clarifying discussions. 



[1] 

[2] 
[3] 
[4] 



[5 



[9 
[10 

[11 
[12: 

[13 

[14 

[is; 

[16: 

[17 
[18 

[19 



[20 

[21 
[22 
[23 

[24 
[25 

[26 
[27 

[28 
[29 



184526 (2009). 
P. Fulde, Phys. 
Eremin, G. Zwicknagi, 
Rev. Lett. 101, 187001 



Y. Kamihara et al., J. Am. Chem. Soc. 130, 3296-3297 
(2008). 

Ren Zhi-An et al.. Chin. Phys. Lett. 25, 2215 (2008). 
Cao Wanget al., Europhys. Lett. 83, 67006 (2008). 
A.V. Chubukov, D. Pines, and J. Schmalian, A Spin 
Fluctuation Model for d-Wave Superconductivity; D. 
Manske, I. Eremin, and K.H. Bennemann, Electronic 
Tfieory for Superconductivity in high-Tc Cuprates and 
Sr2Ru04, K.H. Bennemann and J.B. Ketterson Editors, 
Volume XL Superconductivity: Novel Superconductors, 
Springer- Verlag Berlin Heidelberg (2008). 
P. Thalmeier et al.. Superconductivity in Heavy Fermion 
Compounds, A. Narlikar (Ed.) Frontiers in Superconduct- 
ing Materials, Springer Verlag, Berlin, (2005). 
W.K. Park and L.H. Greene, J. Phys.: Condens. Matter 
21 103203 (2009). 

R.S. Gonnelli et al., Phys. Rev. B 79, 
J. Chang, I. Eremin, P. Thalmeier, 
Rev. B 75, 024503 (2007); L 
P. Thalmeier, P. Fulde, Phys. 

(2008) . 

I.L Mazin and J. Schmahan, Physica C 469, 614 (2009). 
V. Cvetkovic and Z. Tesanovic, Europhys. Lett. 85, 
37002, (2009); 

Physica C 469, (2009), Special Issue on Pnictides. 

L. Benfatto, M. Capone, S. Caprara, C. Castellani, C. 

DiCastro, Phys. Rev. B 78, 140502(R) (2008). 

E. Z. Kuchinskii, M. V. Sadovsku, JETP Letters, 89, 156 

(2009) . 

L. Benfatto, E. Cappelluti and C. Castellani, Phys. Rev. 
B 80, 214522 (2009). 

G.A. Ummarino, M. Tortello, D. Daghero, R.S. Gonnelli, 
Phys. Rev. B 80, 172503 (2009). 

G.A. Ummarino, M. Tortello, D. Daghero, R.S. Gonnelli, 

to be published in J. Supercond. Nov. Magn. 

D.S. Inosov et al.. Nature Physics 6, 178 (2010). 

Ar. Abanov, A. V. Chubukov, and J. Schmahan, Journal 

of Electron Spectroscopy and Related Phenomena 117- 

118, 129 (2001). 

A. Akbari, I. Eremin, P. Thalmeier, and P. Fulde, 
Phys. Rev. B 80, 100504 (2009) M.M. Korshunov and 
I. Eremin, Phys. Rev. B 78, 140509 (2008); T.A. Maier 
and D.J. Scalapino, Phys. Rev. B 78, 020514 (2008). 
I.I. Mazin, D.J. Singh, M.D. Johannes, M.H. Du, Phys. 
Rev. Lett. 101, 057003 (2008). 

P. Popovich et al., Phys. Rev. Lett, 105, 027003 (2010). 

G. M. EUashberg, Sov. Phys. JETP 11, 696 (1960). 

L. Boeri, M. Calandra, I.I. Mazin, O.V. Dolgov, F. Mauri, 
Phys. Rev. B 82, 020506 (2010). 

H. Ding et al., Europhys. Lett. 83, 47001 (2008). 

D. Daghero et al., Phys. Rev. B 80, 060502(R) (2009); 
R.S. Gonnein et al., Physica C 469, 512 (2009). 
M. Tortello et al., Phys. Rev. Lett, 105, 237002 (2010). 
J. Paglione and R.L. Greene, Nature Physics 6, 645 

(2010) . 

I. I. Mazin, private communication. 

G.A.Ummarino and R.S. Gonnelh, Physica C 328, 189 
(1999). 



* Electronic address: 
I mail:giovanni.ummarino@infm.polito.it| 



